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A bottom  reflectivity  measurements  program  has  been  undertaken  to 
determine  the  acoustic  reflection  coefficient  for  specific  areas  of 
ocean  bottom  as  a function  of  incident  angle  and  frequency.  Pulsed 
single  frequency  signals  are  used  for  these  measurements.  Deep  sub- 
merged source  and  receivers  are  oriented  so  that  the  direct  and  first" 
order  bottom  reflected  pulses  are  separable  in  time  and  are  not  interfered 
with  by  the  first-order  surface  reflected  pulses.  The  travel  time  for 
the  propagation  of  both  the  direct  and  bottom  reflected  pulses  is  accurately 
determined  and  recorded  along  with  the  desired  acoustic  signal.  Since 
multiple  receiving  hydrophones  are  used,  the  data  is  also  recorded  in 
a multiplexed  form.  The  requirement  for  efficient  and  expedient  data 
analysis  along  with  the  complexity  of  the  measurement  has  necessitat  d 
the  writing  of  a digital  computer  program.  This  program  computes  the 
reflection  coefficient  and  angle  of  incidence  .for  each  reception  of  data. 

The  program  corrects  for  difference  in  path  length, taking  into  account 
the  depth  of  the  source,  receivers, and  water,  the  bottom  slope,  the 
velocity  gradient  and  the  recorded  travel  time  for  each  reception. 

This  memorandum,  therefore,  discusses  (l)  the  acquisition  and  processing 
of  the  acoustical  data;  (2)  the  definitions  of  the  reflection  coefficients; 
(3)  the  theoretical  relationships  that  properly  determine  the  angle 
of  reflection  and  reflection  coefficient,  which  are  governed  by  the 
geometric  configuration  of  the  source  and  receivers  and  (4)  the  computer 
program  along  with  the  input  and  output  requirements. 
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ACQUISITION  AND  PROCESSING  OF  THE  ACOUSTIC  DATA 

It  is  the  purpose  of  these  studies  to  obtain  a measure  of  the  loss 
that  an  acoustic  signal  undergoes  upon  reflection  from  the  ocean  bottom 
as  a function  of  incident  angle  and  frequency.  This  can  be  accomplished 
by  resolving  the  multipath  structure  of  Figure  1,  and  comparing  the 
amplitudes  of  the  direct  and  the  first-order  bottom  reflected  pulses, 
corrected  in  terms  of  propagation  loss  for  the  difference  in  path  length. 


Velocity 
Prof i le 


Figure  1 Simplified  Ray  Paths 


The  acquisition  of  the  bottom  reflectivity  data  is  planned  around 
the  availability  of  an  AGOR-type  research  vessel  which  has  the  capability 
of  lowering  an  acoustic  projector  to  depths  in  excess  of  10,000  feet. 

This  capability  necessitated  the  development  of  a deep  submergence  package 
which  contains  two  free- flooding  magnetostrictive  scrolls  used  as  projectors. 
An  accurate  depth  sensing  device  is  also  contained  on  the  package.  It 
is,  therefore,  possible  to  lower  the  projector  below  the  sound  channel 
axis.  With  the  receivers  positioned  below  the  sound  channel  axis,  the 
angle  of  reflection  can  be  calculated  with  a greater  degree  of  accuracy, 
and  a greater  signaL-to-noise  ratio  can  be  expected  for  the  bottom 
reflected  pulses  over  the  entire  range  of  incident  angles. 
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One  of  the  desired  parameters  in  this  study  is  the  angle  at  which 
acoustic  energy  is  reflected  from  the  ocean  bottom.  To  obtain  this 
parameter  with  a sufficient  degree  of  accuracy  it  is  necessary  to  measure, 
as  closely  as  possible,  the  travel  time  of  a transmitted  pulse  from 
the  projector  to  the  hydrophone  over  both  the  direct  and  first-order 
bottom  reflected  paths. 

In  order  to  obtain  the  required  measure  of  travel  time,  two  synchronized 
Time/Frequency  Standards  are  used.  A master  is  installed  aboard  the 
transmitting  vessel  and  a slave  at  the  receiving  station. 


Figure  2 Synchronizing  System 


The  master  controls  the  transmission  of  the  projector  while  the  slave 
indicates  the  exact  time  of  the  transmission  at  the  receiving  station 
and  triggers  the  timing  and  control  circuits  to  process  the  data. 


The  desired  data  are  composed  of  pulsed  single  frequency  signals  that 
are  received  at  the  hydrophone,  amnlified,  full  wave  rectified  and 
detected.  The  resultant  signal,  the  envelope  of  the  received  pulse, 
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is  converted  into  digital  form  at  a rate  of  1,000  samples  per  second. 
When  considering  the  processing  of  this  type  of  data  via  a digital 
computer,  it  is  mandatory  that  only  the  necessary  signals  for  proper 
analysis  be  fed  into  the  computer.  It  is  wasteful  both  in  time  and 
money  to  process  signals  or  noise  that  is  not  pertinent  to  the  analysis. 
It  is  necessary,  therefore,  that  the  processing  system  be  accurately 
synchronized  such  that  analog  signals  are  converted  to  digital  form 
only  when  the  desired  data  are  being  received;  thus  a "signal  aperture" 
is  established. 

When  the  slave  Time /Frequency  Standard  indicates  the  exact  time  of 
transmission,  pulses  also  generated  by  the  Standard  are  accumulated 
in  a preset  counter  until  a predetermined  preset  number  is  reached. 

The  choice  of  the  preset  number  is  made  by  an  operator  who  bases  his 
setting  of  the  counter  dependent  on  the  position  of  the  previous 
pulse  in  the  signal  aperture.  When  the  preset  number  is  reached,  a 
pulse  is  generated  by  the  preset  counter  which  is  used  to  trigger  an 
oscilloscope  that  has  its  sweep  rate  set  at  a value  such  that  the 
received  pulse  is  encompassed  on  the  oscilloscope  face.  At  the  exact 
time  that  the  oscilloscope  is  triggered,  the  Franklin  analog-to-digital 
converter  commences  processing  of  the  received  signals.  The  Franklin 
converts  the  signals  for  a predetermined  period  of  time,  dependent 
on  the  pulse  length  of  the  received  pulse.  Thus  the  pulse  is  physically 
observed  in  the  signal  aperture  as  it  is  simultaneously  being  converted 
to  digital  form.  It  is  necessary,  however,  to  be  properly  timed  and 
to  have  all  equipment  functions  properly  sequenced  before  actually 
recording  or  converting  the  desired  signals  to  digital  form. 


Since  the  preset  counter  starts  accumulating  pulses  from  the  Time/ 
Frequency  Standard  at  the  exact  time  of  transmission,  and  controls  the 
conversion  of  the  received  pulses  when  the  preset  number  is  reached, 
the  value  of  the  preset  number  is  proportional  to  the  travel  time  up 
to  the  signal  aperture. 


T = N 

r ~r 


(1) 


T^  - Travel  time 

N - Number  set  into  counter 

f = Frequency  of  pulses  being  fed  into  counter 

The  numerals  set  into  the  preset  counter  are  recorded  through  a coding 
matrix  i nto  the  ODD  records  of  the  digital  magnetic  tape  format  used 
for  these  measurements. 
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Figure  3 Digital  Tape  Format 


The  EVEN  records  contain  the  digital  equivalent  of  the  received  analog 
pulses.  The  travel  time  in  the  ODD  records  is  in  milliseconds  and  is 
the  measure  of  the  time  to  the  beginning  of  the  signal  aperture.  To 
determine  the  exact  travel  time  to  the  signal  in  the  aperture  it  is 
necessary  to  count  the  number  of  digital  samples  taken  up  to  the 
beginning  of  the  received  signal.  The  beginning  of  the  received 
signal  is  defined  as  the  point  at  which  a predetermined  threshold  level 
is  exceeded.  The  number  of  digital  samnles  is  corrected  to  its  equiv- 
alent time  in  milliseconds  and  is  added  to  the  travel  time  number  in 
the  ODD  record.  Thus,  the  exact  travel  time  is  preserved  along  with 
the  data  for  each  recention  of  data. 

The  first  ODD  record  (I.D.  Record)  of  the  taoe  format  contains  ten 
numbers  that  are  coded  for  pertinent  information  about  the  file  of  data. 
It  is  Ttsed  as  identification  of  the  data,  and  to  cross  reference  input 
cards  containing  parameters  used  in  processing  the  data  through  the 

romrvter. 
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The  acquisition  equipment  can  be  operated  in  two  modes  dependent  on  the 
separation  in  time  between  the  direct  and  reflected  pulses.  If  the 
pulses  are  separated  by  100  milliseconds,  two  separate  apertures  will  be 
used  to  obtain  the  data  in  the  EVEN  record,  one  for  the  direct  pulse 
and  one  for  the  reflected  pulse.  These  apertures  are  controlled  by  the 
accumulators  of  a dual  preset  counter.  Thus  the  time  of  conversion  is 
kept  to  a minimum  and  only  the  pertinent  data  is  converted.  When  this 
mode  is  used,  there  are  two  sets  of  5 binary  coded  decimal  numbers 
recorded  in  the  ODD  record;  the  first  is  the  travel  time  to  the  direct 
aperture  and  the  second  is  the  travel  time  to  the  bottom  reflected 
aperture.  If  the  separation  in  pulses  is  under  100  milliseconds, 
then  only  one  aperture  is  used,  and  this  aperture  contains  both  the 
direct  and  reflected  pulses.  Therefore, there  will  be  only  one  set  of 
5 binary  coded  decimal  numbers  recorded  in  the  ODD  record  indicating 
the  travel  time  to  the  signal  aperture. 

The  Franklin  analog-to-digital  converter  system  can  multiplex  input 
signals.  With  the  sample  rate  used  for  these  measurements  it  is  possible 
to  have  up  to  five  channels  multiplexed  in  the  EVEN  records.  It  is 
also  an  option  of  the  system  to  record  data  from  one  hydrophone  through 
two  channels,  one  high  gain  and  one  low  gain  (i.e.  a 10  db  difference  in 
gain  between  channels).  Thus,  if  the  high  gain  channel  overloads,  the 
low  gain  channel  will  contain  the  undistorted  data.  This  increases  the 
dynamic  range  of  the  system  but  must  be  taken  into  account  when  processing 
the  data  through  the  computer. 

The  operator  of  the  preset  controls  has  the  responsibility  of  establishing 
a threshold  value  that  will  be  used  as  an  input  parameter  to  the  computer 
program.  This  threshold  value  is  constant  for  a group  of  pulses, 
usually  associated  with  a station  at  which  approximately  100  events  are 
recorded,  and  is  based  on  the  observed  noise  level.  Additional  input 
parameters  are  necessary  for  the  computer  program;  however,  these  will 
be  discussed  later  on  in  the  memorandum. 


COEFFICIENTS  OF  REFLECTION 

Bottom  reflectivity  measurements  are  performed  with  the  requirement  of 
verifying  a theoretical  model  (reference  a)  which  considers  the  bottom 
as  a medium  made  un  of  a number  of  plane  oara 1 le i layers  of  varying 
acoustic  nronerti es . Absorption  of  both  longitudinal  and  shear  waves 
has  been  considered  in  order  to  describe  the  complex  reflection  and 
transmission  prorcrties  of  t.he  ocean  bottom.  It  ’.■/as  necessary,  therefore, 
to  develop  a computer  program  that  would  process  and  present  the  experi- 
mental data  ir.-  a form  suitable  for  comparison  wi  th  the  predictions  of  the 
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model . 

To  describe  the  extent  to  which  the  amplitude  is  reduced  upon  reflection 
from  the  ocean  bottom  it  was  decided  to  have  the  computer  program 
evaluate  three  different  quantities  for  determing  this  reduction  in 
amplitude.  Considering  a pulse  made  up  of  a single-frequency  sine  wave 
which  has  been  allowed  to  reach  a steady  state  before  it  is  terminated, 
and  is  rectified  and  detected,  it  is  possible  to  characterize  it  by 

(l)  the  time  average  of  the  envelope  of  the  signal 


(2)  the  mean-square  value  of  the  envelope  of  the  signal 

Jrj  fSM 

where  p = envelope  function  of  the  pressure 
T = duration  of  the  pulse 

and  (3)  by  the  peak  pressure  which  is  the  maximum  amplitude  of  the 
steady  state  signal.  By  computing  each  of  these  quantities  for  the 
pulse  of  the  direct  arrival  and  the  pulse  of  the  first-order  bottom 
reflected  arrival  and  finding  their  ratio,  three  separate  coefficients 
are  obtained. 

Of  the  three  reflection  coefficients  which  this  program  provides,  the 
one  based  on  peak  amplitude  is  currently  of  greatest  interest,  since 
it  can  be  interpreted  as  the  ratio  of  steady-state  amplitude  of  a 
continuous  signal  before  and  after  reflection,  and  can  thus  be  compared 
with  the  predictions  of  an  existing  model  (reference  a).  This  model 
is  based  on  a steady- state  solution  for  a single  frequency  plane  wave 
incident  on  a multi-layered  bottom.  Thus,  a comparison  using  the  measured 
peak  coefficient  can  be  made  as  long  as  the  bottom  structure  used  in 
the  model  takes  into  account  the  effective  penetration  of  the  pulse  up 
to  the  time  at  which  the  maximum  occurs.  That  is,  the  maximum  amplitude, 
occuring  "t"  milliseconds  from  the  beginning  of  the  pulse,  can  be 
regarded  as  the  maximum  amplitude  of  a continuous  wave  containing 
components  reflected  at  only  those  interfaces  which  the  pulse  has 
traversed  in  "t/2"  milliseconds. 

At  present,  USL  Computer  Program  No.  0134  based  on  this  model  does  not 
provide  theoretical  values  for  the  other  two  coefficients.  However, 
it  is  felt  that  these  two  quantities  are  both  legitimate  measures  of  the 
ocean  bottom's  reflective  prooerties.  The  time  integral  of  the  envelope 
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of  the  signal  divided  by  the  duration  of  the  pulse  gives  us  the  average 
amplitude  of  the  pulse.  The  time  integral  of‘  the  square  of  the  envelope 
divided  by  the  duration  of  the  pulse  is  proportional  to  the  energy  flux, 
that  is,  the  energy  per  unit  area  per  unit  time.  Thus,  the  ratio  between 
the  bottom  reflected  and  incident  values  of  each  of  these  quantities 
indicates  a reduction  in  the  average  acoustic  pressure  for  one  case, 
and  in  the  flow  of  energy  for  the  other  case. 


THEORETICAL  RELATIONSHIPS 

Unequal  path  lengths  for  the  direct  and  bottom  reflected  acoustic 
pulses  required  that  corrections  be  applied  to  the  ratios  so  that  they 
represent  the  ratios  which  would  actually  be  measured  at  the  water- 
sediment  boundary.  Furthermore,  the  angle  of  incidence  at  the  bottom 
was  required  for  each  transmission, so  that  a method  of  computing  this 
angle  with  reasonable  accuracy  had  to  be  devised.  The  following  assumptions 
and  derivations  served  as  the  basis  for  the  mathematical  relationships 
used  in  the  computer  program. 

ANGIE 

The  chief  assumption  here  is  that  a constant  velocity  gradient  exists 
in  every  part  of  the  medium  which  is  traversed  by  the  rays  we  are 
considering.  This  assumption  is  reasonable,  provided  the  source  and 
receiver  are  both  located  below  the  axis  of  the  permanent  sound  channel. 
Figure  4 represents  the  geometry  of  the  ray  paths  under  such  conditions. 

(See  ton  of  Page  9) 

The  approach  used  here  is  similar  to  that  in  reference  (b),  except  that 
we  require  an  explicit  solution  of  the  angle  o ^ without  knowledge  of  the 
source  angle.  From  Snell's  Law,  it  can  be  shown  that  the  path  of  an 
acoustic  ray  in  a medium  having  a single  constant  gradient  is  the  arc 
of  a circle  having  the  form 

(*-xa)V  (rib)*  - (2) 

where  V,.  = sound  velocity  at  the  source, 

g --  the  velocity  gradient,  a constant, 

»s  - the  source  angle  of  the  ray 
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Figure  4 Geometry  of  Ray  Paths 


Specifically,  choosing  the  coordinate  system  as  in  Figure  4,  we  have 

(x-  ( 

for  path  (a),  and 
for  path  (b). 

Equation  (3)  can  be  expressed  in  terms  of  instead  of  0S  by  applying 
Snell’s  Law: 

Vs  A&C  &s  - \JV  J&C.  o(  . ( 

Thus,  equation  (3)  can  be  written 

( X - x„)'  + - *,)*=  -pr  • < 

d 

It  is  now  nossible  to  obtain  an  expression  for  the  angle  from 
equations  (4)  and  (6)  in  terns  of  known  quantities. 
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Differentiating  ( 4)  and  (6)  we  get 

dl  __  __  X - X. 

^ ' r* 


At  x=xB  , y=yB 


_ y — *0  ^ 

jx  " v* 

, these  are  related  to  tan  c/  as  follows : 
X^-X*  I 


xe-x,  __  4. 

- — — — — • — Wn«K 

and 

- W*  , 

5 

where  the  sign  of  the  tangent  is  determined  by  the  quadrant  in  which 
the  ray  lies.  Utilizing  the  fact  that  yQ  and  y-j  represent  the  depth 
at  which  the  velocity  is  zero,  we  have 

q = '"'b  ~ 0 

* %- 

and 

C,  = *zs.  , 

* 36- a. 


from  which  we  obtain 


3b"3-=  3b-3-=  £ - 

Substituting  this  in  equations  (9)  and  ( 10)  and  eliminating  X , we 
obtain  ® 


Xo  — *,  +- 


■ \x2sl''  £*  , 


i . 

We  now  need  two  more  equations  in  order  to  solve  for  tan  0/  . Consider 
equation  (6);  when  x =■  o,  y _ ys.  Thus  X 


v>-_£  - 

4 4;  ‘1  ^ 
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* 


I 


therefore , 


and  thus  equation  (15)  becomes 


r 


v + rv8j  . t.  _ vj 


3‘ 

When  x = R,  y = y ; therefore,  equation  (4)  becomes 

B .x 


(i i-x.r  ■*-  -§ 

Solving  for  (H  - x ),  we  obtain 


■x  j 


A2C  <A 


R-X,  = ± \f£ (ictf  . 


Blit 


thus 


q — Mg  O • 
i >"*  ’ 


fs*-*T= 


V« 


and  equation  (20)  becomes 

x,  = R _ 


Combining  equations  (14)  and  (23)5  we  obtain 

Xo-Rt\ 

Equating  (24)  to  (18),  we  get 


Xo  — R ta/STS  - )£  + a , 

V <y-  r 3 


i/|  *•*-  = R r + 

We  can  ignore  the  si., ns  of  tne  radicals,  since  these  quantities  will 
all  be  squared  in  the  process  of  solving  for  tan  ^ . 


(17) 

(18) 

(19) 

(20) 
(21) 
(22) 
(23) 
(24' 
(25' 


The  result,  of  squaring  equation  (25)  tw  ce  and  combining  terms  (an 
operation  which  is  a little  too  lengthy  to  be  presented  here)  is  the  cubic 


A iu  .^X 


B *f  C h*~'.  J.  -j-  0 - O 


(26) 
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A=  -KVetR 

8 = XH-Z?  *‘-2M+v£jJ 
C = *W[2'£-W'+'£4ffl] 

0=  -Rv,uv^+2tf  RYv4j+Vk')-  - 

Thus,  if  the  gradient,  range  and  velocities  at  the  source,  receiver, 
and  bottom  are  known,  the  reflection  angle  can  be  calculated. 

RANGE 

Although  the  velocities  and  deaths  can  be  measured  with  reasonable  accuracy, 
horizontal  range  from  source  to  receiver  is  difficult  to  obtain.  Uncertainty 
in  the  horizontal  position  of  a deeply  submerged  package  relative  to  the  ship, 
due  to  the  catenary  which  the  cable  assumes,  increases  the  error  introduced 
in  shin  positioning.  Since  acoustic  travel  time  can  be  measured  fairly 
Precisely,  it  was  felt  that  range  computed  from  travel  time,  average  velocity 
and  depth  would  be  sufficiently  accurate. 
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Consider  Figure  5*  Assuming  a single,  constant  gradient  below  the  axis 
of  the  permanent  sound  channel  (the  point  y = o),  the  average  velocity  is 
readily  computed  from  the  relation 

v = (^tVu,;=  (27) 

Thus,  from  the  straight  line  geometry  of  Figure  5> 

R=  , (28) 

Making  use  of  the  travel  time  t over  the  direct  path  h , which  is 
directly  measured  for  each  transmission,  and  the  average  velocity 
over  h as  computed  from  equation  (27),  we  can  determine  the  range 

from  

R-  )-VjV  - ha'  3s)1  r (29) 


PATH  LENGTH 

Since  the  results  of  this  program  are  to  be  compared  with  the  predictions 
of  a model  which  does  not  consider  loss  due  to  propagation  in  the  water, 
corrections  for  the  path  lengths  of  the  direct  and  bottom  reflected 
arrivals  must  be  applied  to  the  measured  reflection  coefficients.  In 
general,  the  path  length  of  a ray  in  a medium  having  a constant  velocity 
gradient  is  given  by  (Reference  c) 


where 

n = an  integer  which  indicates  the  layer  in  which  the 
velocity  gradient  is  constant, 
gn  = the  value  of  the  gradient  in  the  layer, 

- the  velocity  at  the  upper  boundary 

7 

n+1  - the  velocity  at  the  lower  boundary 
V 

x = the  velocity  at  which  the  ray  vertexes. 

The  inverse  cosines  in  the  brackets  must  be  expressed  in  degrees, 
vertex  velocity  can  be  determined  if  the  velocity  and  angle  of 
inclination  of  any  point  on  the  ray  is  Known.  From  Snell's  Law, 

V - 

C Cs--3’L>  * 

Thus,  V>:  for  both  sections  of  the  bottom  path  can  he  calculated,  since 


The 

(31) 
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the  velocity  at  the  bottom  i6  known,  and  the  angle  which  the  ray  makes 
with  the  horizontal,  which  is  just  the  grazing  angle  ©<  , can  be 
computed.  However,  since  no  angle  is  known  on  the  direct  path,  the 
vertex  velocity  is  unknown,  and  this  path  length  can  only  be  approximated 
by  the  straight  line  h in  Figure  5 and  must  be  computed  from 


i\  = vt 


where 


is  the  travel  time  of  the  direct  arrival. 


PROPAGATION  LOSS  CORRECTION 

As  was  mentioned  earlier,  the  reflected  and  direct  amplitudes  must  be 
corrected  for  differences  in  path  length  in  order  that  their  ratio 
accurately  represents  the  reflection  coefficient.  The  appropriate 
corrections  are  made  to  the  ratios  themselves  by  application  of  a multi- 
plicative factor  derived  from  the  relation  given  in  reference  (d) 

(*q.  m _ a.  (S^-S') 


where 


^ --  (%)  fttr) 


1^  = intensity  at  the  beginning  of  the  path, 

* 2 = intensity  at  the  end  of  the  path, 

Sx  = distance  from  the  source  to  the  point  at  which  1^  is 
measured , 

= distance  from  the  source  to  the  point  at  which  I is 
measured,  ^ 

n = an  integer  to  be  determined  by  the  type  of  spreading  involved, 
a = the  attenuation  coefficient  , 

This  is  the  anti  log  of  the  familiar  propagation  loss  equation.  For  our 
purposes,  spherical  spreading  (n-  -2)  and  an  attenuation  coefficient 
given  by 

a~  0,01  f* 


were  assumed. 


Since  the  reflection  coefficients  are  ratios  of  acoustic  pressures,  we 


J 
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must  take  the  square  root  of  both  sides  of  the  above  equation,  which 
gives  us  the  ratio  of  acoustic  pressures  measured  at  two  points  along 


the  path.  Thus,  if  a and  b are  two  points  on  a ray  path,  then 

S,  lrp o-ff+WV) 

Pw~  V 


(35) 


i Image  Source 

Figure  6 

Identification  of  Pressures  and  Paths 
Used  in  Propagation  Loss  Equations 
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Using  equation  (35),  we  obtain  the  following  relations  between  the 
pressure  amplitudes  at  the  various  points: 

p ~ s 10  ' , 

1 A 4 • 


JU_ 

R. 


S.  _ 


(37-a) 


(37 -b) 


^,0-°.Oo0^VS.) 


(37-c) 


Combining  (37-a)  and  (37-b),  we  get 


Pz 


p: 


ft  ~ c,ct(2  * 

Substituting  P in  (37-c)  into  (38),  we  obtain 

d - 'Sx.  h- 

17  C-C*  Py  ‘ 

Thus,  by  computing  the  ratio  C /C  C from  equation  (37),  in  which 

3 1 2 


(38) 

(39) 


frequency  has  units  of  kc/sec  and  the  path  lengths  S. f expressed  in  kyd 
are  determined  from  equations  ( 30 ) and  (32),  the  corrected  reflection 
can  be  obtained.  It  should  be  noted  that  S , which  corresponds  to  the 
distance  between  the  source  and  the  point  at  which  the  acoustic  pressure 
"at  the  source0  is  measured  is  taken  to  be  one  yard  by  convention, 
to  prevent  indeterminate  terms  from  arising. 


COMPUTER  PROGRAM  DESCRIPTION 

USL  Computer  Program  No.  0289  processes  observed  data  for  the  determination 
of  reflection  coefficients. 

The  input  is  in  the  form  of  digitized  magnetic  tape  (Datrac),  and 
Hollerith  cards  as  detailed  in  Table  I.  Two  cards  are  required  for 
each  file  of  information  on  Datrac  Tape.  The  tape  consists  of  up  to 
three  files,  the  even  records  of  which  contain  voltage  readings  to  three 
significant  figures  oroncrtional  to  the  acoustic  pressure.  The  first 
odd  record  in  each  file  is  an  identification  word  which  is  comoared  with 
the  corresuonding  identification  word  from  cards.  The  first  even  record 


16 
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in  each  file  is  calibration  data  and  is  used  to  obtain  a multiplicative 
constant  by  which  low  gain  data  are  modified  to  correspond  to  high  gain 
information,  when  applicable.  Subsequent  odd  records  contain  information 
as  to  the  time  (in  milliseconds)  until  the  aperture  x’or  the  direct  pulse 
information  begins,  and  possibly  the  time  to  the  start  of  the  bottom 
reception  aperture.  If  only  one  value  aopears  in  the  odd  record,  it  is 
assumed  that  the  signal  was  recorded  continuously,  after  the  beginning 
of  .he  direct  aperture,  through  the  bottom  reception.  The  even  records, 
after  the  first  for  each  file,  contain  reception  data  for  direct  end  bottom 
reflection.  The  urogram  allows  for  up  to  five  channels  of  information 
from  these  even  records.  Each  channel  may  pertain  to  one  hydrophone, 
or  two  channels  may  be  allocated  to  one  hydrophone  so  that  high  and  low 
gain  data  are  present.  If  the  high  and  low  gain  condition  exists,  the 
proerarn  has  an  alternative  means  of  getting  results  when  the  maximum 
reception  voltage  is  observed  in  the  high  gain  channel.  It  is  assumed 
that  calibration  data  will  be  present  in  the  second  record  of  each  file 
even  if  the  high  and  low  gain  option  is  not  being  used. 

The  program  has  two  main  divisions.  The  first  pertains  to  obtaining  and 
processing  all  the  necessary  data  from  tape,  and  the  second  deals  with 
the  computation  and  printing  of  the  desired  output  results.  The 
separation  of  these  two  major  parts  is  at  statement  number  200  in  the 
program. 

In  the  second  part  of  the  program  is  a subroutine  for  solving  a specific 
third  degree  polynomial . The  polynomial  is  in  terms  of  tangent  oC 
whereof  is  the  angle  of  reflection  that  the  acoustic  ray  makes  with  the 
bottom.  The  subroutine  uses  the  Newton-Raphson  method  for  obtaining 
one  root  of  the  cubic  equation.  The  other  roots  are  obtained  by 
reducing  the  oolynomial  one  degree  and  solving  the  resulting  equation 
by  the  quadratic  formula.  The  first  root  obtained  is  accepted  if  the 
remaining  roots  are  imaginary,  or  if  the  remaining  roots  lay  outside 
the  range  of  +2C  to  If  one  of  these  conditions  is  not  met,  a 

message  will  be  printed  by  the  computer.  A starting  value  of  alpha 
is  part  of  the  input  data  for  the  program  arid  is  necessary  for  the 
operation  of  this  iterative  process.  Experience  with  the  data  determining 
the  coefficients  of  these  polynomials  has  indicated  that  a choice  of  angles 
between  +60°  ana  +70°  will  insure  convergence  if  the  data  affecting 
direct  travel  time  are  correctly  given  to  the  program. 

In  the  first  part  of  the  program  semi-fixed  criteria  are  used  to 
determine  when  the  direct  and  bottom  pulses  have  begun  and  ended.  In 
the  case  of  the  direct  pulse,  ten  consecutive  samples  above  threshold 
are  sought  to  determine  the  start  of  the  pulse,  and  five  consecutive 
samples  below  the  threshold  voltage  indicate  the  end  of  the  Pulse. 
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The  bottom  reception  starts  when  ten  consecutive  samples  appear  above 
the  threshold  voltage.  The  end  occurs  when  one  sample  falls  below  or 
equals  the  threshold  voltage.  However,  if  subsequently,  five  consecutive 
pulses  exceed  the  threshold  voltage,  then  all  the  samples  which  occurred 
since  the  previous  termination  are  considered  as  part  of  the  bottom 
reception.  In  addition,  the  bottom  reception  includes  all  those  that 
continue  to  stay  above  threshold  for  a predetermined  time  or  until 
the  samdes  again  equal  or  fall  below  threshold,  whichever  comes  first. 

The  Fortran  Program  is  available  as  Appendix  A.  Also,  samples  of  the 
printed  output  are  contained  in  Figures  7 and  8.  Figure  7 represents  the 
output  when  sense  switch  one  is  down,  and  Figure  8 covers  the  print-out 
with  sense  switch  two  down.  These  options  may  be  exercised  singly  or 
jointly.  With  sense  switch  two  down,  all  the  calculated  data  are  put  on 
a master  tape  togeter  with  the  pulse  values  above  threshold  and  all 
input  information.  This  can  then  be  used  for  statistical  studies  done 
by  other  specially  written  programs. 


SUMMARY 

• 

This  memorandum  discusses  an  integrated  system  for  the  acquisition, 
processing  and  analysis  of  Bottom  Reflectivity  Data.  A technique  has 
has  been  successfully  adapted  and  favorable  results  are  presently 
being  obtained.  The  acquisition  of  data  is  semi-automatic,  once 
sychronization  is  established,  and  the  data  is  obtained  in  digital  form 
for  computer  processing.  The  computer  program  also  discussed  has 
been  written  based  on  an  analysis  procedure  outlined  in  this  memorandum. 
This  program  determines  reflection  coefficients  as  a function  of  angle 
of  incidence.  The  analysis  procedure  and  the  output  requirements  of 
the  program  are  such  that  an  effective  comparison  of  the  results  can 
be  made  with  predicted  theoretical  results.  These  results  are  obtained 
from  a multi-layered  theoretical  model  of  the  bottom  via  a computer 
program  where  the  input  parameters  are  determined  from  analyses  of 

s sediment  cores  obtained  at  the  same  area  of  ocean  bottom  where  the 

acoustical  measurements  were  made. 


! 

. 
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35-37  XX. X Approximate  Angle  of  Incidence  in  Degrees 


TABLE  I CONTINUED 


21-30  XXXXXXX.XXX  Digital  Sampling  Rate  ( Sample s/Msc ) 


PROGRAM  OUTPUT  SENSE  SWITCH  ONE  DOWN 


IBM  TAPE  NO.  205 


FILE  NO.  1 


ID  WORD  APPROX.  RANGE 
1060354011  1.0  KIDS. 


EVENT  FREQ. 
A l.OKC 


PUU3E  WIDTH 
10.MSC. 


THRES.  VOLT. 
0.015  VOLTS 


SOURCE  DEPTH 
2925. FT. 


REC.  DEPTHS 
3862. FT. 
4165. 


IBM  TAPE  NO.  285 


FI  IE  NO. 


RECEPTION  NO. 


CHANNEL  NO.  2 


PEAK  COMHJTATIONS 
DIRECT  PEAK 
1.064 


REFIECTED  PEAK 
0.091 


DIRECT  PRESS. /SOURCE  PRESS.  REFLECTED/DIRECT 

0.75780631E-03  0.110 


AVERAGE  COMPUTATIONS 

DIRECT  SUM  DIRECT  AVERAGE  BOTTOM  SUM  1 BOTTOM  SUM  TOTAL  BOTTOM  AVERAGE  BOTTOM  A V. /DIRECT  AV. 

11.680  0.433  1.459  1.459  0.049  0.145 

ENERGY  FLUX  COMPUTATIONS 


MEAN  SQUARED  DIRECT  SIM  OF  BOTTOM  1 SQUARED  MEAN  SQUARED  BOTTOM  BOTTOM/DIRECT 


0.334 

0.094 

0.003 

0.017 

MISCELLANEOUS 

90-ALPHA 

48-90DEG. 

BOT.  S AMPULS 
30 

DIR.  SAMPLES 

27 

FIRST  BOT.  SAMPIES 

30 

DIB.  TRAVEL  BOT.  TRAVEL 

0.810  SEC.  1.045  SEC. 

COMP.  RANGE 
1.280  KYD. 

IBM  TAPE 

NO.  285 

FILE  NO.  1 

RECEPTION  NO.  1 

CHANNEL  NO.  3 

PEAK  COMPUTATIONS 

DIRECT  PEAK  REFLECTED  PEAK  DIRECT  PRESS . /SOURCE  PRESS .-  REFLECTED /DIRECT 

O.896  0.092  0.74205969E-03  0.125 


AVERAGE  COMPUTATIONS 

DIRECT  SUM  DIRECT  AVERAGE  BOTTOM  SUM  1 BOTTOM  SUM  TOTAL 
10.437  0.417  1.741  1.741 


BOTTOM  AVERAGE  BOTTOM  AV. /DIRECT  AV. 

0.050  0.145 


ENERGY  FLUX  COMPUTATIONS 
MEAN  SQUARED  DIRECT 
0.274 


SUM  OF  BOTTOM  1 SQUARED  MEAN  SQUARED  BOTTOM 

0.113  0.003 


BOTTOM/DIRECT 

0.024 


KISCELIANEOUS 

90-ALPHA 

BOT  SAMPLES 

DIR.  SAMPLES 

FIRST  BOT.  SAMPLES 

DIR.  TRAVEL 

BOT.  TRAVEL 

51-67DEG. 

35 

25 

35 

0.(27  SEC. 

1 .006  SEC. 

COMP.  RANGE 
1.280  KYD. 
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FIGURE  7 


PROGRAM  OUTPUT  SENSE  SWITCH  TWO  DOWN 


DATRAC  TAPE  MASTFR  TAPF 

NO.  FILE  RECFPT IONS  NO.  TOTAL  RECORDS 

?16  1 15  155  145 


FVFNT  ID  NUMBFR 

A 1060112114 


HYDROPHONE 

NUMBFR 

1 

2 


DEPTH  IN 
FEET 

1640.0 

1752.0 


QO-ALPHA  (DFG) 
MAX.  MIN, 

52.640  52.308 

51.165  52.777 


FIGURE  8 
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APPENDIX  A 


C DETERMINATION  OF  REFLECTION  COEFFICIENTS  FROM  OBSERVED  DATA  0289 

C 

S  ION  OCT  016060606060 
S IBB  OCT  226060606060 
ICC  OCT  236060606060 
IDD  OCT  2A6060606060 
iEE  OCT  256060606060 
I2Z  OCT  716060606060 

DIMENSION  1DI2 ) .API  2 ) »CALN( 5 ) .CALDI 5 ) . WD ( 5 » 100 ) .DR ( 5 ) * YR ( 5 ) .VRI5 ) . 

1 I DUMP (91*  NSD ( 5 ) . NSB ( 5 ) . T TD ( 5 ) . T TB ( 5 ) *LD < 5 ) .LB ( 5 > * 

2 SUMA  ( 5 ) .SUMS  ( 5 ) .WMAXI  5 ) .SlJMB  ( 5 ) »SUMBS(  5 ) »BMAX(  5 ) «SV  IX<  5 ) 

3 * NSB2 ( 5) .SVl (61 . TAPI  180) .TMP ( 160  I * AMNA ( 5 I * AMXA ( 5 I 
READ  14  » 1 RT  »IFT  »MIT 
CLA  I RT 
TZE  *1 

DO  1111  L = 1 * 1 RT 
RTD  1 

1111  CONTINUE 
IOD 

1 READ  INPUT  TAPE  3 . 5 . I BM . I F I LE . I D3 . I D4 .E VENT .PW .FREQ . IGA  I N .MUL T .X . 

1 NREC F . IWPR.APR ,AANG .FAC ,6 . TH . VA ,DA . DB , DS .DR ( 1 ) 

READ  INPUT  TAPE  3.6  . ( DR ( I ) . I =2 . 5 I , PM 
DO  2 1—1.5 
AMNA 1 1 ) =90.0 

2 CAL N I 1 ) =0.0 

3 FORMAT  ( 15HDATRAC  TAPE  NO. , I 5 . 5X8HF I LE  NO. . I 3 »5X 1 5HMASTER  TAPE  NO. 
1 .A5/46HRECORDS  ON  MASTER  TAPE  TO  START  OF  DATRAC  FILE. 16) 

I DUMP ( 6 ) = I DUMP ( 6 ) +2 

5 FORMAT  ( I3.U.2A5.A1.F3.0.F2.1.2I1.F2.0.I3.U.2F3.1.F2.1.F5.A.F3.3 
1.2FA.0.3F5.0) 

6 FORMAT  ( 4F5.0.F10.3 ) 

CLA  G 
TZE  *12 

PUNCH  3. IBM. IFILE.MI T. IRT 
CALL  DATSP  ( 7 . I D ( 1 ) .0 . 1 .0 . 0 . 0 .0 . 0 ) 

DO  666  1=1.2000 
STZ  CALN(l) 

666  CONTINUE 

NH=MULT/ IGAIN 
CAL  ID3 
ARS  6 
SUB  I D ( 1 ) 

TZE  *10 

7 PRINT  8 

8 FORMAT (38H  ID  ERROR.  CHECK  INPUT  CARDS  WITH  TAPE) 

PAUSE  7007 

10  CAL  I DA 
ARS  6 
SUB  10(21 

TZE  *15 
TRA  *7 

12  NRECF=(NRECF*2)+3 
DO  13  J* 1 .NRECF 
RTD  7 
IOD 

I DUMP ( 6 ) ■ I DUMP ( 6 ) + 1 

13  CONTINUE 


A-l 


UST,  TECH  MEMO 


n n n 


i 


i 

i 


s 

I 


\ 


1A  FORMAl  I 2 I 5 * A5  ) 

I DUMP (61  = ! DUMP < fa )-2 

GO  TO  1 

1*<  IF  < SENSE  SWITCH  It  1 6 » 2 1 

16  PRINT  2 0*1 8M  . I F 1 LE . ID3* IDA.APR.EVFNT ,FREO.Pw.TH,DS.(DR<I)»l*l.NHI 
? 0 FORMAT  ! 1 HI  / 1H0  * 1 8X  1 ?H  I BM  TAPE  NO  . . I A . 9X  8HF  I |_E  NO. » I 2 / 1 HO » 1 X9 1 HI  0 

1 WORD  APPROX,  RANGE  EVENT  FRtU.  PULSE  WIDTH  THRES.  VOLT.  SOU 
2RCE  DEPTh  RE'.  PEP T HS/ 1 X . 2 A 5 . 3X ,F 5 . 1 . 5HK Y DS . * 5X . A 1 . 3X , FA. 1 , 2HKC * 

3 T'UF5.0*i<HMSC.  . 3X.F5.3.6H  VOL  IS . 5X »F 6 • 0 • 3HFT . *3X»F7.0*3HFT./82X. 
4F7.0/82X,F7.0/82X.F7.0/82X.F7.0> 

21  TH=10.*TH 
VS=DS-DA 
Y6=DB-DA 
VS=YS*G+VA 
VB=YB*G+VA 
I M=  1 
IE*0 
NCA  = 1 00 

I 5-IWPR-I 100*MULT ) 

TPL  *38 
NCA-IWPR/MULT 
I WA  = 0 
IREM=0 
TRA  *39 

38  IWA=IWPR/( MULT  *100 ) 

I REM* ( IWPR-IWA*MULT*100)/MULT 

39  I SAVE 7 = I W A 

I SAVE8= IREM 
I SAVE9=NCA 


USE  CALIBRATION  DATA 


CALL  DATSP  (7.WDI 1*1) *0*0* NCA  »MUL  T ♦ 5 * I * 0 1 
IF  ( IGAIN-i 1 A6.A6.A0 
AO  DR ( A ) =DR ( 2 ) 

DRI 3) =DR(2  1 
DR  I 2 1 *DR I 1 ) 

DO  A5  1*1 .MULT  * I GA I N 
CALNI I 1=0.0 
CALD( 11=0.0 
X=  I +1 

DO  AA  J = 6 * 1 5 

CALNI I 1 =CALN I 1 l+WDI I.J) 

CALDI I )=CALD( I 1+WDIK.J) 

AA  CONTINUE 

CALNI K)=CALN( 1 1 /CALDI I 1 
AS  CONTINUE 

A6  IF  (SENSE  SWITCH  21  AA6.A7 

AA6  WRITE  OUTPUT  TAPE  1 * AA7 * 1 BM . I F I LE . 1 ON . NH , I D3 . IDA . APR , E VENT .FREQ* 

1 PW.TH.DS. I DR ( I 1 . 1 = 1 .5 1 .CALNI 2 1 .CALM  A) 

IRT=IRT+1 

A A 7 FORMAT  I 13.  I 1.5XA1 , I 1 .2A5.FA. l.AltF3.1»FA.0.FA.2.F5.0»5F6.0»2E15.8 
1 1 


A 7 DO  A8  J* 1 .MULT 
YRI J)=DR( Jl-DA 
A8  VR I J 1 =YR I J ) *G+VA 
50  API  1 1 =0.0 
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non  (/icn</'t/io-.  n nn 





API?) -0.0 
N S B 1 = 0 
l T - 1 
l.HN'-O 

LMN 1 >0 

OBTAIN  APERTURES 


12  = 1 

CALI  DATS P (7.API1)  .0.0. 1.2. 2.1*1) 

1 E * I E + 1 

I DUMP  t 6 ) = I DUMP  ( 6 ) T-2 
I V T = I R T 
63  CIA  API?) 

TIE  *255 
C LA  API  1 ) 

TZE  *65 
TRA  *75 

65  API  1 ) =AP( 2 ) 

AP I 2 ) “0.0 

75  CALL  DATSP  ( 7 , WD (1 . 1 > . 0 .0 , NCA ,MUL T , 5 . I 2 .0 > 


SEARCH  FOR  START  OF  DIRECT  PULSE 


< 

i 


i ! 

: * 


5 

i 

r 


TV=TH 
I 3 I M 
I CD3 1 
NS  = 0 
N TV  = 0 

DO  81  J» 1.100 
NS=NS+1 

IF  (WD(I.J)-TV)  80.00.79 

79  NTV=NTV+1 

IF  (NTV-10)  81.83,33 

80  NTV  = 0 

81  CONTINUE 
JFK  3 JFK  + 1 

IF  (JFK-3)  8000.8000.9000 
8000  BACKSPACE  7 
GO  TO  75 

9000  SUMAI I 1=0.0 

IF  (SENSE  SWITCH  1)  9001,190 

9001  PRINT  82.1 

82  FORMAT  I 1H0 , 1 1HCHANNEL  NO..I3.A7H  SHOWS  NO  SlUNAL  ABOVE  THRESHOLD 
1 I N 100  SAMPLES) 

GO  TO  190 

83  NS=NS-10 
JFK  =0 

T TO  I I ) =NS 

T TO  C I )=TTD<  I )/PM+AP  ( 1 >+X 
I 1 =NS+ 1 
l.0 I I ) =0 

Lam  -o 

84  SUMA ( i ) =0.0 
SUMS ( I » =0.0 
\S1=0 
NTV*0 


A-3 
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NSD (1  t sO 
15  = 0 

SUM6X  =0 • 0 

SUM8SX-0.0 

WMrtX  ( ) ) = WD  I 1 * I 1 ) 

in«o 

START  OF  DIRECT  PULSE 

86  DO  93  J*Il  ,NCA 

IF  ( WD ( I » J )-9«  99  J 07.9A.94 
b 7 IF  I T V-WD I l » J > 1 39.08.88 

88  NTV=NTV+1 

SUM8X  =SUM8X  +WD ( I . J) 

SUM9SX=SUM8SX+(WD( 1 ,J)«*2> 

TMP(LT)=WD( I »JI 
lt=lt  + 1 

IF  (NTV-5)  92.97,9? 

89  IF  <WD<  I ,J  J-WMAXU)  ) 91,91.90 

90  W'MAX  ( 1 ) =WD(  1 .J) 

91  SUMA { I ) =SUMA ( I M-WDl I .JI+SUMBX 

SUMS  < 1 ) “SUMS  < I ) + ( WD ( 1 . J)**2  M-SUMBSX 
CL A SUM3X 
TZE  *992 
DO  991  L=1.NTV 
L TM  = L T-NTV+L— 1 

991  TAP(LTM)=TMP(LTM) 

992  SUMBX c0 • 0 
N T V*Q 

SUMB  SX«0  • 0 
TAP ( LT ) *WD ( I »J) 

IF  (LT-22)  2992. Z663. 2662 

2662  LT 1 “LT-22 
L T =22 

2663  IF  (SENSE  SWITCH  2)  266A.2669 

2 66A  WRITE  OUTPUT  TAPE  1 , 998  . I BM . I F i LE . I E , 1 . I BB , ( T AP ( L ) . L = 1 . LT ) 
LMN=LMN+1 
L T = 1 
CLA  LT  1 
TZE  *92 

DO  2665  L * 1 » LT 1 
2665  TAP(L )=TAP(L+22 » 

LT*LT1+1 
L T 1 = 0 
GO  TO  92 
2669  LT1*0 

lt=o 

2992  LT=LT+1 

92  NSD ( 1 ) *NSD ( I 1 + 1 

93  CONTINUE 
GO  TO  113 

DIRECT  PULSE  SWI ICHED  TO  LOW  GAIN  CHANNEL 

9 A LT  = 1 

CLA  LMN 
TZE  * 99A 

A-U  usl  tech  memo 

9i3_U_6,t? 
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in  in  in  in  t/i  c/i  in  in  in  ooo 


00  993  L 3 1 » LMN 
BACKSPACE  1 

993  CONTINUE 
LMN*0 

994  IF  (IGAIN-1)  197,197,95 

95  IF  (LDIII-1)  995,197,197 

995  TV  3 TV/CALN ( I +1  I 
T TO ( 1+1 ) *TTD  I I ) 

SUMAI  I I *9.99 

LDI I ) *1 
1*1  + 1 
LDI I)*l 

IF  l I CD-I)  84,84,96 

96  BACKSPACE  7 
ICD*1 

12*1 

CALL  DATSP  I7,WD(1,1) ,0,0,NCA, MULT, 5, 12.01 
I 1*NS+1 
GO  TO  84 

DIRECT  PULSE  IS  TERMINATED 

97  I1=NSDI I >-15+1+11 
T V*TH 

15*0 
LT-LT-6 
TZE  *999 

996  IF  (SENSE  SWITCH  2>  997,999 

997  WRITE  OUTPUT  TAPE  1 ,998 , IBM, I F I LE , I E , I , I BB , ( T AP ( L ) »L«1 ,LT ) 
LMN=LMN+1 

998  FORMAT  ( I 3 , 1 1 , I 3 ♦ 1 1 . 1XA1 .22F5.2 ) 

999  LT-1 

I RT» I RT+LMN 

LMN1-LMN 

NSDI I ) *NSD I I > -4 

NV*NSDI I ) 

CLA  LDI I > 

TZE  *98 

SUMAI I ) *SUMA( I )*CALN< I > 

SUMS! I )*SUMS( I ) * I CALN I I >**2> 

WMAXI I l-WMAXI I ) +CALN I I t 

98  CLA  AP I 2 ) 

TZE  *103 

1 2* I WPR/ ( 2*MULT ) + 1- 1 CD* 100 
TMI  *102 

NCA-100 

BACKSPACE  7 

1 1* I WPR/2+1 

I 2* I I WPR/200)*100+1 

I 1*  I 1 1-12) /MULT  + 1 

I W4* I IWPR-I2+1)/IMULT*100) 

I REM* I I I WPR- 1 2+1 ) - 1 MULT* I W4*l 00 ) J/MULT 
CLA  IW4 
TZE  *99 
TRA  *100 

99  NCA* IREM 

100  CALL  DATSP  I 7. WD 1 1 , 1 ) »0 ,0, NCA .MULT , 5 , 1 2 ,0 ) 


A-5 


USE  TECH  MEMO 

9I3-L-65 
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I CD  = l 
GO  TO  104 

102  I 1 = 100+12 
12  = 1 

GO  TO  104 

103  NS  1 *5 

in/,  i 1 1 u 

N T V*0 

IF  (Il-NCA)  105,105,125 

SEARCH  FOR  START  OF  BOTTOM  PULSt 

105  DO  112  J = I 1 * NCA 
NS1 =N5 1+1 

I F (WD(ltJ)-TV)  107,107,106 

106  NTV=NTV*1 

IF  (NTV-10)  112.136,136 

107  NTV*0 

112  CONTINUE 
GO  TO  125 

GET  MORE  DATA  FOR  DIRECT  PULSE 

113  I 2= I 2+ ( 100*MULT ) 

I 5 = NSD( I ) 

IF  (ICD-IW4)  115.114,180 

114  CLA  IREM 
TZE  *180 
NCA* IREM 

115  BACKSPACE  7 

CALL  DATSP  (T.WDIl.l) ,0*0, NCA, MULT, 5, 12,0) 

I CD* I CD+1 
11*1 

GO  TO  86 

GET  MORE  DATA  IN  SEARCH  FOR  BOTTOM  PULSE 

125  IF  (IC.D-IW4)  127,126.184 

126  CLA  IREM 
TZE  *184 
NCA- I REM 

127  I 2* I 2+ ( 100*MULT ) 

I 5-NS1-5 
BACKSPACE  7 

CALL  DATSP  I 7 , WD ( 1 . 1 ) .0 .0 . NCA ,MULT , 5 , I 2 ,0 ) 

II  = 1 

ICD-ICD+1 
GO  TO  105 

128  NCA=100 

START  OF  BOTTOM  PULSE  IS  AT  THE  END  OF  PREVIOUS  DATA 

I 2 = I 2-( 100*MULTI 
I CD* I CD-I 
I 1*100+I 1 
BACKSPACE  7 

CALL  DATSP  ( 7 , WD( 1 • 1 ) ,0 ,0 , NCA ,MULT ,5,12,0) 

A-6  USL  TECH  MEMO 
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GO  TO  140 


START  OK  BOTTOM  RECEPTION 

136  NS  1 *NS 1-10 
CLA  AP  ( 2 ) 

TZE  *138 
T TB ( I ) *NS 1 
GO  TO  139 

138  NTV=AP I 1 >»PM 

T TB I I )*NS1+NS+NV+NTV 

139  T TB ( I » - T TB I I ) /PM+AP ( 2 > +X 
11*11 +NS1-1 5-5 

TMI  *128 

140  SUMBC I )*0.0 
SUMBSI I »“0.0 
SUMBX*0*0 
SUMBSX-0.0 
NS8 ( 11*0 
NSB2I I ) *0 
NSB 1 *0 
NTV*0 
NUM-1 

BMAX ( I >*NV 
NFSD-FAC*BMAX<  I I 
BMAX ( I ) * WO  < 1*11) 

I SAVE1* 1 1 

ISAVE2-I2 

1 SAVE3* ICO 

ISAVE4-NCA 

LT*  1 

LT1*0 

LMN*0 

SVI I 1 )*0.0 
SVIX( I )*0.0 

I 5* ( I WPR- 1 2+1) /MULT -1 1 + 1 
IF  (15  - NFSD)  141*142*142 

141  NFSD* I 5 

142  15*11 +NFSD-1 

IF  (15-100)  145,144.144 

144  15*100 

145  NFSD*NFSD+I 1-1-15 
DO  159  J* 1 1 * 1 5 

IF  (WD( I.J)-9.99)  150*175.175 

150  IF  (TV-WD(I.J))  154.151,151 

151  CLA  NUM 
TZE  *152 

SVI ( I ) *SUMB ( I ) 

SVI X C 1 ) *SUMBS( I ) 

NSB2 ( I ) *NSB ( I ) 

NUM*0 

152  NTV* 1 

153  SUMBX*SUMBX+WD( I *J) 

SUM8SX-SUMBSX+ (WD< I .J)**2) 

NSB 1 *NSB1+ 1 
TMP(LT)*V<D(I.J) 

lt*lt+i 

A -7  USI,  TECH  MEMO 
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GO  TO  159 

lr’A  IF  I WD<  i . J >-8MAX  ( I »)  156,156,155 
155  BMAXI  1 ) *WD I I ♦ J) 

S 156  CLA  N T V 

S TZE  #158 

NTV=NTV+1 

IF  (NTV-5)  153,153.157 

157  N T V = 0 

N SB ( I ) = NSB ( I > +NSB 1 
SUMB ( 1 ) = SUMB ( I ) ♦SUMBX 
SUMBS  t I I *SUMBS( 1 M-SUMBSX 
S CLA  MSB  1 

S TZE  #958 

DO  957  L* 1 »NSB1 
LlM-LT-NSBl+L-1 

957  TAP(LTM)*TMP<LTM) 

NS81*0 

958  SUMBX  *0*0 
SUMBSX=0.0 

158  SUMB I I ) =SUMB<  f )+WD( I ,J) 

SUM8S ( I J«SUMBS( I )♦< WD( I ,J>**2  > 

NSB ( I ) *NSB ( I )+l 
TAP(LT)*WOU»J) 

I BG*  1 
IND*22 
L T 5*0 

2958  IF  (LT-22)  6959.1959.1959 
1959  IF  (SENSE  SWITCH  2)  2959,8959 

2959  WRITE  OUTPUT  TAPE  1 » 998 , I BM» I F I LE , I E , I . I CC , ( T AP ( L ) »L* 1 BG , I ND ) 
LMN=LMN+1 

LT  *LT-22 
L T 5*  1 

I8G*IBG#22 
I ND* I ND+22 
GO  TO  2958 
S6959  CLA  LT5 

S TZE  *8959 

S CLA  LT 

S TZE  #8959 

00  7959  L* 1 » LT 
LT 1 * I ND+L-22 
7959  TAP(L)-TAP(LT1) 

8959  LT*LT+1 

159  CONTINUE 
C 

C GET  MORE  DATA  FOR  BOTTOM  RECEPTION 

C 

S CLA  NFSD 

S TZE  *190 

IF  (ICD-IW4)  162,161.190 
S 161  CLA  IREM 
S TZE  #190 

NCA* I REM 

162  I2*I2+( 100*MULT) 

BACKSPACE  7 

CALL  DATSP  ( 7 , WD (1,1) ,0 »0, NCA .MULT ,5,12,01 
1 CD* I CD# 1 

A -8  USL  TECH  MEMO 
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I la) 

GO  TO  142 


BOTTOM  PULSE  SWITCHED  TO  LOW  GAIN  CHANNEL 


175  L T * 1 
NSB  1 » 0 

CLA  LMN 
TZE  *2175 
DO  1175  L*  1 *LMN 
SAC  K SPACE  1 
1175  CONTINUE 
l MN  = 0 

2175  IF  (IGAIN-ll  197,197.176 

176  IF  < L B < 1 » - 1 I 177,197,197 

177  LB( II -1 

TTBI 1+1 >*TTB( I ) 

TV*  T V/CALN ( 1 + 1 ) 

1 = 1+1 
LB  ( I 1 * 1 

IF  ( I CD-I  ) 140,140.178 

178  BACKSPACE  7 
ICD-ISAVE3 

I 1 = I SAVE  1 
12*1 SAVE2 
NCA-ISAVE4 

CALL  DATSP  ( 7»WD( 1, 1 ) ,0,0. NCA, MULT ,5,12,01 
GO  TO  140 

180  IF  (SENSE  SWITCH  1)  181.187 

181  PRINT  182.1 
IRT  * 1 RT  + LMN 

182  FORMAT  ( 1H0, 24HRECORD  ENDED  FOR  CHANNEL . 14, 2X39HWI TH  DIRECT  PULSE 
1ST  I LL  ABOVE  THRESHOLD) 

GO  TO  187 

184  IF  (SENSE  SWITCH  1)  185,187 

185  PRINT  186,1 

186  FORMAT  < 1HO,24HRECORD  ENDED  FOR  CHANNEL , 1 4 , 2X43HW I TH  NO  BOTTOM  BOU 
1 NCE  PULSE  ABOVE  THRESHOLD) 

187  I * I M 

BMAX ( I )-0.0 
SUMS ( I 1*0.0 
SUMBS ( I )*0.0 
SVIIII *0.0 
SVI X( I 1-0.0 
T TB ( I >-0.0 
NSB ( I )»C 
NSB2 ( I ) -0 
LB( I 1 -0 
LB( I+l)-0 
L T 5 * L MN 
LT-1 
NSB 1 *0 
LMN-0 
CLA  LT5 
TZE  *190 
DO  189  L-1.LT5 
189  BACKSPACE  1 

A-9  USL  TECH  MEMO 
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IRT=JRT-lT5 


END  OF  TAPE  WORK  FOR  ONE  CHANNEL 

190  CLA  LB  < I ) 

TZE  * 1 9 A 

SUMBI I ) *SUM8 < I >*CALN( I ) 

SUMBS i 1 ) *SUM3S( 1 t *(  CALN( I ) ** 2 ) 

SVK1  >*SVl  ( I )*CALN<  I ) 

SV1X(1)*SVIX(I ) * ( CALN ( I >**2) 

BMAXI  l I * 8MAX ( 1 )*CALN(  I ) 

194  IM*IMfiGAlN 
L T = L T- 1-NSB 1 
TZE  *9196 

IF  (SENSE  SWITCH  2)  195.9196 

195  WRITE  OUTPUT  TAPE  1 .998 . I BM. I F I LE . I E . I . I CC , ( T AP < L > . L-l * LT ) 
LMN*LMN+1 

9196  L T=  1 

IRT*IRT>LMN 
NSB  1*0 
LMN*0 
LMN1-0 
12*1 

I W4* I SAVE7 

IREM-ISAVE8 

NCA-1SAVE9 

IF  (IM-MULT)  196.196.201 

196  BACKSPACE  7 
GO  TO  75 

197  IF  (SENSE  SWITCH  1)  198.200 

198  PRINT  199.1 

199  FORMAT  ( 1 HO  * 20HV0LT AGES  FOR  CHANNEL . I A . 2X7 1HHAVE  EQUALED  OR  EXCEED 
1 ED  THE  MAXIMUM  VOLTAGE.  THIS  CHANNEL  WAS  ABORTED.) 

200  SUMA ( IM) *0.0 
CLA  LMN1 

TZE  *194 

DO  9201  L*l.LMNl 
9201  BACKSPACE  1 
IRT* IRT-LMN1 
GO  TO  194 

ALL  TAPE  WORK  COMPLETE.  PREPARE  FOR  FINAL  COMPUTATIONS 

201  1*1 

202  IM* I 
CLA  SUMA ( I ) 

TZE  *245 
J*LD( I ) ♦LB ( 1 ) 

IF  (J-l)  220.203.219 

203  J*m 
CLA  LD( I ) 

TZE  *204 
TRA  *207 

204  NSD ( J ) *NSD( I ) 

TTD ( J ) *TTD( I ) 

SUMA ( J) -SUMA ( I ) 

SUMS ( J ) "SUMS ( I ) 


A 10 


USL  TECH  MEMO 


CO  CO  to  Co 


wMAX ( J 1 =WMAX ( 1 ) 

GO  T G 219 
207  NS»(  J>*NSBf  1 ) 

NSB2 ( J > =NSB2 I 1 ) 

T TB ( JJ-TTBl  l ) 

SUMBl J)=SUMB( 1 ) 

SUMP S ( J ) =SUMBS I 1 1 
BMAX  < J > =BMA  X ( I ) 

SVI  ( J ) *> SV  1 ( I ) 

SV I X ( J ) *SV I X ( I ) 

219  1*1-1 

220  SO*. 001 

RANG* SQRTF II  SO*  T TD  ( I ) * I G/2 . 0 * < YR  ( 1 > 4 yS  H-VA  ) > **2- < V R ( I )-YS >**2 ) 

V*  VR ( 1 ) 

BANG* AANG 

221  CALL  CUBIC  ( BANG. G, RANG. VB.V, VS. ANG) 

CLA  ANG 

TZE  *253 
VX*VB/COSF  < ANG) 

S1*S0*VX/  ( G*  3-.  0 ) 

S2=S1* I ATANF ( SQRTF  < VX**2-V**2 ) /VI-ANG) 

S 1*S 1 * < ATANF ( SQRTF ( VX**2~VS**2 ) /VS ) -ANG) 

S3  = SO*SQRTF(RANG**2-HDR(  1 ) -DS  ) **2  ) / 3 .0 
F3=-.0005*FREQ**2*(S3-S0) 

F3*10.**F3«S0/S3 

RAT I0=.0O05*FREQ**2*( S1+S2-S3) 

RATIO  * ( (S1+S2 )*10.**RAT10)/S3 
RRT  0=»RAT  IO*6MAX  ( I )/WMAX<  I ) 

RBSDS=RATIO**2*SUMBS< 1 )/SUMS( I ) 

SUMA ( I ) *SUMA< I ) / 1 0 • 

AB«NSD I I 1 

AD * SUMA ( 1 ) /AB 

SUMS ( I l •SUMS l l ) / ( 100* *AB ) 

SV! ( I >*SVI ( 1 )/10. 

SUMB ( I ) * SUMB ( I ) / 1 0 • 

WM A X ( I ) »WMAX (11/10. 

BMA  X ( I ) »PMAX ( I )/10. 

TT0( i )«S0*TTD< I ) 

TT8( 1 >-S0*TT6< I ) 

AB  =NSB ( I ) 

SUMBS ( i > *SUMBS ( 1 )/( 100.*AB) 

AB  = SUMB ( I ) /AB 
RBD»RAT10*AB/AD 
DC  T 
NOP 

svm  i i-svixi  i )/ioo. 

ANG=90.-ANG*57.2957795 
RANG«S0*RANG/3.0 
IF  (SFNSF  SW I T rw  ii 

222  PRINT  225. IBM, IFILE.IE.I 

225  FORMAT  ( 1HO/1H0.3X12HIBM  TAPE  NO. . I A, 7X8HF I LE  NO. , I 2 . 7X 1 3HRECEPT 10 
IN  NO. » I A » 7X 1 1HCHANNEL  NO. .12) 

PRINT  226  »WMAX ( I ) ,3MAX( I) .F3.RRTD 

226  FORMAT  ( i.H0,i7HPEAK  COMPUT  AT  IONS/5X,  1 1HDIRECT  PEAK  . 7X 1 AHREFLECTED 
1PEAK .6X27HDIRECT  PR  ESS. /SOURCE  PRESS . » 6X 16HREF LECT ED/U I RECT /6X ,F6 . 
23.1 5X.F 5.3*1 7X.E 15.8. 16X.F6. 3 ) 

PRINT  227 .SUMA ( I ) .AD.SVl ( I ) .SUMBl 1 ) .AB.RBD 

A -11  USL  TECH  MEMO 

91 3 -4-65 


CO  </>  CO  CO 


227  FORMAT  < 1H0.20HAVER AGF  C.OMPUT  A T I ON  .*  / 5x  10HD 1 REC  T SUM , 3X 14HD I RECT  AV 

1ERAGE  .5X1/HBOTTOM  1 . 5 X 1 6HBOT  TOM  SUM  TOTAL  .5X1  4HB0T  T OH  AVERAGE* 

25X2 1H6CT TOM  AV. /DIRECT  AV. /6X  »F8  • 3 . 10X  . F 5 . 3 . 12X  ,F  6 . 3 . 1 1 X »F8.  3 * 13X  » 
3F6.3»16X,r7»3) 

PRINT  2 2 8*  SUMS l I ) *SV I X ( l ) .SUMBSI  1 ) .RBSDS 

22 8 EORMAI  I 1HU.2AHENERGY  FLUX  COMPUT AT  I ONS/5X2 1H  MEAN  SQUARED  DIRECT 
1.5X23HSUN  OF  BOTTOM  1 SQUARED  * 3X2 1H  MEAN  SQUARtD  BOTTOM  .3X13HBOTT 
20M/DIRECT/11X.F8.3.19X.F8.3.19X.F8.3.15X.F7.3) 

PRINT  c29  .ANG.NSB  l I ) .NSDl  i ) »NSB2  m.TTD(l).TTBU)  • RANG 

229  FORMAT  I 1HG.13HMISCELLANEOUS/5X8H90-ALPHA.5X12HBOT.  SAMPLES » 5X 12HD 
UR.  SAMPLES*  5X13MFJRST  bOT  . SAMPLES . 5X l 1HDI R . TRAVEL  * 5X I1HBOT . TRA 
2VEL  .5X11HCOMP.  RANG  E / AX  . F 6 . 2 . 4HDEG.  , 7X  , I 4 , 1 3X  * I A * 1 6X  . I 4 , 1 3X  . F 7 . 3 . 
34HSEC. » 5X  *F  7 » 3 * 4HSEC  < , 5X . F 7 . 3 . 4HKYD. ) 

230  IF  (SENSE  SWITCH  2)  231.245 

231  CLA  BMAX ( I ) 

TZE  *245 

WRITE  OUTPUT  TAPE  1 . 2 32  . I BM , I F I LE . I F , I , 1 DD , WMAX ( 1 ) .BMAX ( I ) . RRTD . 

1 SUMA ( 1 ) .AU.SVI ( I ) « SUMB ( I ) .AB.RBD.NSBI  I ) .NSD<  I ) . F 3 

WRITE  OUTPUT  TAPE  1 . 2 32  . I BM , I F I LE . I E . I . I EE . SUMS ( I ) »SVl X ( I J , SUMBS ( I 
1 ) .RBSDS.TTDI  H .TTB( I ) . RANG . ANG . DUMMY ,NSB2 ( I ) 
lRT=IRT+2 

2 32  FORMAT  ( I 3 , I 1 . I 3 , I 1 . 1 X . A 1 , 9F9 . 4 . 2 1 4 , E 1 5 . 8 ) 

M*0 

233  IF  (IGAIN-T)  240.240.234 

234  IF  (1-2)  235.236.237 

235  AMXAI  1 + 1 ) = AMXA ( I ) 

AMNA ( 1+1 ) * AMNA ( I ) 

GO  TO  240 

236  AMXAI l-l )*AMXA< I ) 

AMNA ( 1-1 ) a AMNA ( I ) 

GO  TO  240 

237  I F ( I - 3 ) 2 38.235.236 

238  STOP  1010 
240  CLA  M 

TNZ  *245 

AMXAI  1 t .MAX  IF ( AMXAI I)  .ANG) 

AMNA I l ) *M1 N1F ( AMNA ( I ) .ANG) 

M=  1 

GO  TO  233 

245  I*IM+IGAIN 

IF  I I -MULT ) 202.202.246 

246  IM=>1 

IF  (SENSE  SWITCH  1)  250.9997 

250  PRINT  251 

251  FORMAT  ( 1H1  ) 

9997  IF  (SENSE  SWITCH  5)  9998.9999 

9998  I DUMP  < 1 ) * — 3 
IDUMP ( 2 > * + 0 
IDUMP(3)=2*IFILE 
IDUMP ( 4 ) «-7 
IDUMP ( 5 ) =+0 
IDUMP ( 7 > *-l 

1 DUMP ( 8 ) * + 0 
I DUMP ( 9 ) ■ I R T 
CALL  DUMP ( IDUMP ) 

9999  CONTINUE 

IF  ( I R T-20000 I 9095.9090.9090 

A 12  USL  TECH  MEMO 
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9090  PRINT  9091. 1RT 

9091  FORMAT  ( 1H 1/  1H0 .2  5H  i APE  UNIT  1 IS  FULL  WITH  .15.85H  RECORDS.  CORE 
1WILL  BE  DUMPtD  ON  TAPE  UNIT  5.  RtSlART  WITH  A NEW  TAPE  ON  TAPE  UNI 
2T  1./1H1) 

WRITE  OUTPUT  TAPE  1.276.IZZ 
ENO  FILE  1 
REWIND  ] 

!RT*0 

GO  TO  9998 

9098  IF  ( I E-NRECF ) 50.270*270 

252  FORMAT  l 1H0. 1AHAPPR0X.  ANGLE*  .F6.2 . 10X 13HRECEPT ION  NO. .15) 

253  PRINT  252. BANG. IE 
BANG* bANG ♦5.0 

IF  (BANG-AANG-5.0)  221.221.25* 

25*  IF  (IVT-IRT)  285.2*5.2*5 

255  CALL  DATSP  <7.WD(l.l)»0.0.NC A, MULT *5*12.0) 

IF  ( SENSE  SWITCH  1)  256.2*6 

256  PRINT  259.1 

259  FORMAT  (1H0.35HBOTH  APERTURES  FOR  RECEPTION  NUMBER . I *. 1 XA5HARE  ZER 
10  AND  THIS  RECEPTION  HAS  BEEN  SKIPPED.) 

GO  TO  2*6 

270  RTD  7 
100 

I DUMP ( 6 ) * I DUMP ( 6 ) ♦ 1 
IF  (SENSE  SWITCH  2)  271.1 

271  PRINT  280.IBM. IFILE.IE.MIT.IRT.EVENT.ID3.ID* 

00  1271  L*1.NH 

m«l+igain-i 

PRINT  3271*L.DR(M). AMX A ( M ) « AMNA ( M ) 

1271  CONTINUE 

DO  1272  L*1 . 5 
AMXA ( L ) *0  * 0 

1272  AMNA ( L ) *90 .0 

IF  (IFILE-IFT)  1.272.1 

272  PRINT  275.IRT. IFILE 

275  FORMAT  ( 1H1/1 HO »1*H THERE  ARE  NOW  .I5.*0H  RECORDS  ON  TAPE  UNIT  1 AF 
ITER  COMPLETING. I3.22H  FILES  ON  TAPE  UNIT  7./1H1) 

PUNCH  1*  • I RT 

WRITE  OUTPUT  TAPE  1.276.IZZ 

276  FORMAT  (9XA1) 

ENO  FILE  1 
GO  TO  1 

280  FORMAT  ( 1H0/1H0/1H0/1H0.9X11HDATRAC  TAPE . 18x1 3HMASTER  TAPE/5X3HN 
10..2X16HFILE  RECEPT IONS. 12X3HNO. .AX13HT0TAL  RECORDS/ *X I3.5XI1. 
26X13. 1*XA5»8XI5/1H0/3X5HEV ENT . 5X9H ID  NUMBER / 5XA 1 , 7X . 2A5 / 1H0/ 5X2 1HH 
2 YDROPHONE  DEPTH  I N»7X1*H90-ALPHA  ( DEG ) /7X6HNUM3ER . 7XAHFEET ♦ 

*8X*HMAX. . 8X*HMI N* ) 

285  MA* I VT+1 

DO  287  L-MA.1RT 
BACKSPACE  1 
287  CONTINUE 
IRT-IVT 
IM*1 

PRINT  290 

290  FORMAT  ( 1H0/1H0.26HCOMPLETE  RECEPTION  SKIPPED) 

GO  TO  250 

3271  FORMAT  (10XI1.8XF6.1 . 6XF6. 3 . 6XF6. 3 ) 

END (0.1. 1.1.1) 
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